home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / cungtr.z / cungtr
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCUUUUNNNNGGGGTTTTRRRR((((3333FFFF))))                                                          CCCCUUUUNNNNGGGGTTTTRRRR((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CUNGTR - generate a complex unitary matrix Q which is defined as the
  10.      product of n-1 elementary reflectors of order N, as returned by CHETRD
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CUNGTR( UPLO, N, A, LDA, TAU, WORK, LWORK, INFO )
  14.  
  15.          CHARACTER      UPLO
  16.  
  17.          INTEGER        INFO, LDA, LWORK, N
  18.  
  19.          COMPLEX        A( LDA, * ), TAU( * ), WORK( LWORK )
  20.  
  21. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  22.      CUNGTR generates a complex unitary matrix Q which is defined as the
  23.      product of n-1 elementary reflectors of order N, as returned by CHETRD:
  24.  
  25.      if UPLO = 'U', Q = H(n-1) . . . H(2) H(1),
  26.  
  27.      if UPLO = 'L', Q = H(1) H(2) . . . H(n-1).
  28.  
  29.  
  30. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  31.      UPLO    (input) CHARACTER*1
  32.              = 'U': Upper triangle of A contains elementary reflectors from
  33.              CHETRD; = 'L': Lower triangle of A contains elementary reflectors
  34.              from CHETRD.
  35.  
  36.      N       (input) INTEGER
  37.              The order of the matrix Q. N >= 0.
  38.  
  39.      A       (input/output) COMPLEX array, dimension (LDA,N)
  40.              On entry, the vectors which define the elementary reflectors, as
  41.              returned by CHETRD.  On exit, the N-by-N unitary matrix Q.
  42.  
  43.      LDA     (input) INTEGER
  44.              The leading dimension of the array A. LDA >= N.
  45.  
  46.      TAU     (input) COMPLEX array, dimension (N-1)
  47.              TAU(i) must contain the scalar factor of the elementary reflector
  48.              H(i), as returned by CHETRD.
  49.  
  50.      WORK    (workspace/output) COMPLEX array, dimension (LWORK)
  51.              On exit, if INFO = 0, WORK(1) returns the optimal LWORK.
  52.  
  53.      LWORK   (input) INTEGER
  54.              The dimension of the array WORK. LWORK >= N-1.  For optimum
  55.              performance LWORK >= (N-1)*NB, where NB is the optimal blocksize.
  56.  
  57.      INFO    (output) INTEGER
  58.              = 0:  successful exit
  59.              < 0:  if INFO = -i, the i-th argument had an illegal value
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.